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We discuss recent work on the transverse momentum dependent sea quark density and its 
application to forward Drell-Yan production. 



1 Introduction 

Scattering processes with a single hard scale are well described in QCD within the framework 
of collinear factorization. The treatment of multi-scale processes, on the other hand, is more 
involved. In this case, generalized factorization formulas are needed^ to gain control over large 
logarithms in higher orders of perturbation theory. Such formulas typically involve transverse- 
momentum dependent (TMD), or "unintegrated", parton distribution and parton decay func- 
tions. A broad class of multiple-scale events is given by small- x processes. These are one of the 
main sources of final states in the central region at the LHC, and lead to sizeable rates of forward 
jet production at the LHC0H At small x, TMD parton distributions arise naturally as a conse- 
quence of high energy factorization and BFKL evolution 4 . /cT-factorization - - provides then the 
matching of these high energy factorized TMD distributions to collinear factorized distributions. 
For Monte Carlo applications a convenient description is given in terms of the CCFM evolution 
equation^ which interpolates for inclusive observables between DGLAP and BFKL evolution^. 
This therefore supplies a natural basis for a Monte-Carlo realization of fc^-factorization, such as 
that provided by the Monte Carlo event generator Cascade^. 

Computational tools based on TMD parton densities have so far been developed within a 
quenched approximation where only gluons and valence quarks are taken into accountEEl While 
this captures correctly the leading contributions at small x, it is mandatory to go beyond this 
approximation in order to include preasymptotic effects and to treat final states associated with 
quark-initiated processes such as Drell-Yan production. 

In this contribution we present work 11 in this direction, and its application to forward 
Drell-Yan production. For further detail we refer toEH 



2 Definition of a TMD sea quark distribution and off-shell qq* — > Z coefficient 



The unintegrated sea-quark distribution is analyzed in^ to logarithmic accuracy a s (a s In x) n 
based on the off-shell TMD gluon-to-quark splitting function^. This is obtained by generalizing 
the expansion in two-particle irreducible kernels of^l to finite transverse momenta, and reads 
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Here A = q — zk with k and q transverse momenta of the off-shell gluon and quark respectively, 
while z is the fraction of the 'minus' light cone momentum of the gluon which is carried on by 
the t-channel quark. Although evaluated off-shell, the splitting probability is universal. Once 
combined with the gluon Green's function, it takes into account the small x enhanced transverse 
momentum dependence to all orders in the strong coupling. In this approach the transverse 
momentum of the sea quark arises as a consequence of subsequent branchings at small x, with 
no strong ordering in their transverse momenta. 




(a) (b) 

Figure 1: (a): If the vector boson is produced in the forward region, the sea quark density becomes 
sensitive to multiple small x enhanced gluon emissions, leading to a fc^-dependent gluon density (b): 
Schematic factorization of the partonic qg* —> Zq process of a) into the g* — > q* splitting and the 

qq* — S- Z coefficient. 

To relate this parton splitting kernel to forward vector boson production, we analyze the 
flavor exchange process g*q — >■ Zq, see Fig. [TJ At high (partonic) center of mass energy, this 
process can be treated according to the "reggeized quark" calculusl 14 l 15 [ The latter extends the 
effective action formalism^, currently explored at NLO^, to amplitudes with quark exchange 
in terms of effective degrees of freedom, the so-called reggeized quarks [HLL21 The use of the 
effective verticest 14 ! 15 ! ensures gauge invariance of the coefficients relevant to perform the high- 
energy factorizationlMl for vector boson production, despite the off-shell parton. 

If taken literally, the reggeized quark calculus leads for the g*q —> Zq process to a rather 
crude approximation to the g* —> q* splitting function, associated with the lightcone momentum 
ordering condition which sets the 'plus' momenta of the off-shell quark for the g* —> q* splitting 
to zero. For Eq. ([l]) this corresponds to the limit z — >• 0. It is however possible^ to relax this 
kinematic restriction and to keep z finite, while maintaining the gauge invariance properties of 
the original vertex. For the g* — > q* splitting this yields then precisely the splitting function 
Eq. 0. 

On the other hand, in the qq* — > Z coefficient the high energy limit sets the 'minus' compo- 
nent of the quark momentum to zero. It proves to be possible to relax the ordering prescription 
also in this case. It is thus interesting to investigate the effect of these kinematic corrections, 
which are subleading in the collinear and high energy limits. In^ we express the off-shell 



coefficient for the Z-boson cross section as 



ff M ._>z = V2G F M 2 z (V q 2 + A 2 q )—S(zxix 2 s + T- M 2 Z ). (2) 

Here the variable T parametrizes the off-shellness of the i-channel quark. In the collinear limit 
T -)■ so that Eq. Q agrees with the lowest order qq — > Z coefficient. For the general off- 
shell case, T interpolates between the squared transverse momentum of the off-shell quark, if 
strong minus momentum ordering is fulfilled, and modulus of the four-momentum transfer, if 
this condition is relaxed. Correspondingly, the qg* —> qZ cross section is expressed in terms of 
convolutions in transverse momentum and four momentum transfer^ respectively. 

3 Numerical analysis 

Fig. [2] shows a numerical comparison 11 of the factorized formulas discussed above with the 
qg* — > qZ matrix element result and with an expression which uses only the collinear splitting 
function. For small |A|, the differences between t and /cr-factorized expressions are numerically 




Figure 2: (a): A 2 dependence of the differential cross section da/dA 2 for small |A|: (solid) full; (dashed) no 
plus-momentum ordering; (dot-dashed) no plus-momentum and minus-momentum ordering; (dotted) collinear 
approximation. All but the last curve overlap in this region. We set X1X2S = 2.5A/§, k 2 = 2 GeV 2 . (b): Relative 
deviations in the differential cross section da/dA 2 : (dashed) no plus-momentum ordering; (dot-dashed) no plus- 
momentum and minus-momentum ordering. 

small, and both expressions are close to the full result; as |A increases, we find that the 
deviations due to the kinematic contributions by which the two expressions differ become non- 
negligible, and that the i-factorized expression gives a better approximation to the full result. 

Future extensions of the above results concern large- x contributions l^ 21 ! 22 ! 23 ! 24 -!^ parton 
shower Monte Carlo implementations^ inclusion of full quark emissions in the evolution, see^l 
for related work in the context of next-to-leading order BFKL evolution. 
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